Experimental quantum Zeno effect in NMR with entanglement-based measurement 
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We experimentally demonstrate a new dynamic fashion of quantum Zeno effect in nuclear magnetic 
resonance systems. The frequent measurements are implemented through quantum entanglement 
between the target qubit(s) and the measuring qubit, which dynamically results from the unitary 
evolution of duration Tm due to dispersive-coupling. Experimental results testify the presence of "the 
critical measurement time effect", that is, the quantum Zeno effect does not occur when Tm takes 
the some critical values, even if the measurements are frequent enough. Moreover, we provide a first 
experimental demonstration of an entanglement preservation mechanism based on such dynamic 
quantum Zeno effect. 

PACS numbers: 03.65.Ta, 03.65.Xp, 76.60.-k 
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I. INTRODUCTION 

The quantum Zeno effect (QZE) describes the situation 
of the inhibition of transitions between quantum states 
by frequent measurements [1]. As observed in some ex- 
periments, e.g., using trapped ion [2], cold atoms [3], cav- 
ity quantum electrodynamics [4], nuclear magnetic reso- 
nance (NMR) [5], QZE was often regarded as the ex- 
perimental witness of projection measurement, or called 
the wave-packet collapse (WPC), since the first interpre- 
tation of QZE was made according to WPC. However, 
many people questioned this opinion by re-explaining 
these experimental observations with some dynamic fash- 
ions without invoking WPC [6]. In the approach of 
Ref. [7], each measurement is implemented as a dynami- 
cal unitary evolution driven by a dispersive-coupling of 
the measured system to the apparatus with duration t,„. 
Actually, the free evolution causes the deviation of the 
system from its initial state, while the dynamic measure- 
ments can interrupt the evolution by adding a phase fac- 
tor to the resulting state of the system, leading to QZE. 
However, the dynamical phase effect depends on the mea- 
surement time Tm- When Tm takes some critical values 
Tm, each dynamical phase factor in the measurements 
corresponding to some integer phase in units of 27r, the 
system is unaffected by the measurements. In this case, 
no QZE occurs. We call this phenomenon the critical 
measurement time effect [7]. 

In this paper, we experimentally reveal this r™- 
dependence in a NMR ensemble when the measurements 
are treated by unitary dynamical processes. We first car- 
ried out experiments with one single-qubit system and 
one measuring qubit. Besides the effects predicted by 
the conventional QZE, the role of critical is also 



clearly demonstrated in the experiments. From the view 
of the experiment, the dynamic measurement model is 
more compatible with the physical reality in compari- 
son with the projection measurement in respect of the 
QZE. Therefore it can be regarded as an active mecha- 
nism protecting the system from deviating from its initial 
state. We also experimentally implement such a scheme 
of entangled-state-preservation in a two-qubit system via 
QZE, which is significant to quantum information and 
computation. 



II. A DYNAMIC APPROACH FOR QZE AND 
THE CRITICAL MEASUREMENT TIME 
EFFECT 

A general dynamic approach for QZE is described 
with a sequence of N frequent measurements inserted 
in a unitary free evolution described by U (r), as illus- 
trated in Fig. 1(a). Each measurement is of duration 
Tm with equal time intervals t + Tm, and is also de- 
scribed by a unitary operator M = M {Tm) = e 
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The corresponding Hamiltonian 
Hm = Hs + Ha + Hint describes the time evolution 
of a closed system formed by the measured system S 
plus the apparatus A with the Hamiltonian Hs and Ha, 
correspondingly. M dynamically results in an entangle- 
ment between S and A with the initial state \a). Here, 
\sj) are the orthonormal states of S, but the states \aj) 
of A need not to be orthonormal with each other. The 
entanglement-based measurement means that one could 
readout the system state \sj) from the apparatus state 
\aj). Such measurement M is realized by the disper- 
sive coupling Hint oi A to S, which is so strong that 
we do not consider the relatively weak free evolution 
during measurements. Usually we chose \sj) to be the 
eigenstate of the system Hamiltonian Hs, which satisfies 
[Hs,Hint] — 0, thus M obviously represents a quantum 
nondemolition (QND) measurement [8]. 
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FIG. 1: (color online), (a) General dynamic approach for 
QZE and (b) the corresponding schematic diagram for ex- 
periments. Frequent measurements M are driven by the in- 
teraction Hamiltonian Hm, while the free evolution U of the 
system are implemented by RF pulses applied on the system. 
In the end, we read out the information about the system 
just from the apparatus via QND measurement. The free- 
evolution [/-process is neglected during strong measurement 
M-process. 



As illustrated in Fig. 1(a), the total time evolution de- 
scribing the QZE reads as Utot{t) = [M (r™) f/ (r)]^ 
with the total free evolution time t = Nt fixed. In the 
limit T and meanwhile t keeps constant, Utot {t) 
becomes diagonal with respect to \sj), thus such evolu- 
tion process due to dispersive-coupling inhibits the tran- 
sitions among the states \sj). Remarkably, Utot{t) is 
not always diagonal for arbitrary Tm even in the limit 
r — > 0. When is accessing certain critical values 
T*^ ~ 27m/A, n = 0,1,2,..., the system exhibits no 
QZE, where A roughly represents the energy level spac- 
ing of the system, and thus the precise form of t*^ de- 
pends on the concrete model. 

Figure 1(b) shows the corresponding schematic di- 
agram for the experiment to demonstrate the above 
entanglement-based QZE in a NMR spin system. The 
system consists of nuclear spins = 1,2, ...n), while 
the apparatus is nuclear spin Jq, called the measuring 
qubit. In a static longitudinal magnetic field, the (n+l)- 
spin system has the natural Hamiltonian 

n n 
3=0 j<k,=0 

where Vj is Larmor frequency of spin j and Jjk is the 
scalar coupling strength between spins j and k. The 
free evolution of the system is implemented by radio 
frequency (RF) pulses, and in multiply rotating frame 
U (t) = Q-^Hfr^eT jg then driven by the Hamiltonian 



J = l 



(2) 



Here the chemical shifts Sj = Vj — and we assume 
the individual nuclear spin Ij can be independently ex- 



cited with frequency fJj (e.g. hereto-nuclear NMR sys- 
tems) and the strength Pj of the RF field is so large that 
the spin-spin couplings of strength Jij could be ignored 
when the RF fields are applied to the system (usually 
the hard RF pulses have this property Pj ^ \Jij\)- A 
sequence of RF pulses with frequencies flj are period- 
ically applied to the system. Between the RF pulses, 
the spin-spin couplings are employed to implement an 
entanglement-based measurement A/ ~ AI {Tm) through 
a measurement Hamiltonian Hm'- 



H 



M 



oj 



(3) 



(i.e., So = 0), and 



= 277 J2 JojIzH- Due to the fact Pj > |Jy|, the 



Here, we have chosen Ha 

n 

H ^ " 

role of U (r) in the total time evolution is dominated with 
respect to M in our experiment, which is just opposite to 
the scheme illustrated in Fig. 1(a). However, this does 
not influence our result as we only require the U and 
M processes can be well separated in time domain. In 
order to observe the dynamic QZE with the QND mea- 
surement, we require to prepare the initial state |s) of 
the system into an eigenstatc of Hs- 

To ascertain whether the system deviates from its 
initial state at the end of time evolution, we need to 
read out the information about the system through a 
QND measurement on the measuring qubit. The in- 
teraction Hamiltonian Hint in Hm is chosen to satisfy 
[Hint,Hs] = 0, thus M represents a QND measure- 
ment. Furthermore, [Ha^ Hint] = guarantees a mea- 
surement scheme analogy to the Ramsey interference, 
which can be used in our experiments to detect the de- 
viation of the final state from its initial state. To this 
end, the apparatus Iq is prepared in a superposition state 
I a) = (|0) + |1)) /V2. Then Utot (t) will evolve the ini- 
tial product state |v(0)) ~ \s) §5 \a) to a quantum en- 
tanglement \ip{t)) = -j^{U- \s) |0) + U+ \s) |1)), where 

(t) = {Optot (t) |0) and U+ (t) = {l\Utot (*) |1). U± {t) 
has the similar limitation behavior limjv->oo \U± (t)\ 1 
as Utot (t)-, thus the system is frozen in its initial state by 
the frequent measurements. Consequently, we can ob- 
tain the state information of the system by measuring 
the magnitude of the conherence of spin Iq, i.e., the off- 
diagonal element of its reduced density matrix: 



V 



{s\Ulit)U^it)\s) 



(4) 



When the QZE occurs, U± (t) freezes the initial state 
|s) up to a change of the phase factor, thus V equals to 
unity. However, when U± (t) evolves the system away 
from its initial state, then V should present an oscillat- 
ing dynamics. Accordingly, the behavior of V provides us 
the information of whether the QZE happens. Addition- 
ally, the experimental values of T> can relatively easily be 
obtained, thanks to the quadrature detecting technology 
in NMR signal detection. 
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III. EXPERIMENT 

The experiment was realized at room temperature on a 
Bruker AV- 400 spectrometer {Bq = 9.4T) . The physical 
system we used is the ^^C-labeled Diethyl- fluoromalonate 
dissolved in ^i^-labeled chloroform [10], where ^^C is cho- 
sen as the measuring qubit /q, and ^^F and as the sys- 
tem qubits Ii and I2, respectively. The J-coupling con- 
stants J12 = 48.3 Hz, J02 = 160.7 Hz and Joi = -194.4 
Hz. We first initialized the system in a pseudopure state 
(PPS) pooo = (1 - e)l/8 + e |000) (000| using the spatial 
average technique [11], where 1 representing the 8x8 
identity operator and e w 10~^ the polarization. Then 
the measuring qubit (^^C) was prepared into the super- 
position state I a) = (|0) + |1)) /a/2 by a tt/2 pulse along 
y axis. In order to observe the dynamic QZE in differ- 
ent quantum systems, i.e, a single-qubit system and a 
composite two-qubit system, we performed two sets of 
the experiments. Figure. 2 shows the experimental pulse 
sequences for observing the entanglement-based QZE in 
the singlc-qubit system and the composite two-qubit sys- 
tem. 
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FIG. 2: (color online). Pulse sequences for observing 
entanglement-based QZE: in a single-qubit system (a), and 
in a composite system for a product state |00) (b) and for a 
pseudo-entangled state (c). Each M-process lasts about 

2.5 ~ 3.0ms for the single-qubit system and about 28 ~ 33ms 
for the composite system. Here, ^=5° ~ 6°, 77 = 4007rTm, 



di 



2IJ12I 



, d2 



2| ■^01 



and ds 



2IJ02I 



To observe QZE in a single-qubit system, we took 
only spin Ii as the system while decoupling spin I2 
during the whole experiment. The initial state of spin 
Ii, \s) = |0), and the natural Hamiltonian Hnmr hi 
Eq. (1) can be the measurement Hamiltonian Hm, i-C, 
Hm = 2ty5iII + 2ttJqiI^II. We can obtain the critical 
measurement time = n/ (5i + Joimo) (n = 0, 1, 2, . . . ) 
for this case. In the experiment, we set 5i ~ SOOiJz, 
the strength of the RF field Pi = 18000i7z > Joi and 
the duration of these RF pulses r = 1/is. A series of 
RF pulses and measurements were performed on the sys- 
tem for the different measurement times. We measured 
the NMR signal intensity of ^^C as a function of the 
pulse number N , which is directly proportional to the 
coherence T) of spin /g. Figure 3(a) shows the exper- 
imental data for the measurement time at the critical 



value Tm = T*j = 2.517tos (denoted by the dots), at 
= 2.55ms nearby the critical value (denoted by the 
triangles) and at t„i = 3ms far from the critical value 
(denoted by the squares), along with the theoretical ex- 
pectations (denoted by the solid lines) obtained by nu- 
merical simulations. As expected, when the measurement 
time Tjn is set at the critical value rj^, NMR signal inten- 
sity presents a Rabi oscillation; when Tm is far from rj^, 
the usual QZE occurs where the frequent measurements 
would largely inhibit the unstable system from evolving 
to other states; when t,„ is close to r^, we observe the 
deviation from usual QZE and T) shows a suppressed os- 
cillation. This r,„-dependent entanglement-based QZE 
now is shown in our NMR experiment. 

Meanwhile, we observed the decay of the NMR sig- 
nals in Fig. 3(b) (left side). This is mainly caused from 
the relaxation and the RF inhomogeneity. The trans- 
verse relaxation time are respectively about 300ms for 
^^C and 800ms for ^^F, while each M-process lasts about 
2.4 ~ 3ms (Fig. 2). In order to improve experimental pre- 
cision and get a better observation of our phenomenon, 
we engineered the unitary operation [A/(t,„)[/(t)]'^ of 
the k cycles as a single shaped pulse by the gradient as- 
cent pulse engineering (GRAPE) algorithm [12], where k 
ranges from 1 to N . All these GRAPE pulses are of the 
duration around 2.5ms with the theoretical fidelity above 
0.99. Thus the decay caused by relaxation can be almost 
neglected. We can see the deviation is maximize when 
reaches to rj^, ranging almost from to 1. These pulses 
are also designed to be robust against the RF inhomo- 
geneity. The GRAPE-bascd results shown in Fig. 3(b) 
(right side) give good description of the entanglement- 
measurement QZE. 

In the second set of the experiment, we further demon- 
strate entanglement-based QZE in a composite system 
consisting of spins Ii and I2 ■ Here, we considered two dif- 
ferent initial states of the system: the product PPS |00) 
and the pseudo-entangled state |^+) = (|01) -I- |10)) /^/2 
[13] obtained from |00) by NOT gate, Hadamard gate 
and CNOT gate. Therefore, two kinds of spin-spin 
couplings are involved in the measurement Hamiltoni- 

ans for QND measurements. For the product state, 

2 2 
Hm = 2ttY^ 5jll + 2tt Jjklil^ (i-e., the natural 

j = l j<k,=0 

Hamiltonian of these three-spin NMR system in the mul- 
tiply rotating frame). In this case, the critical measure- 
ment time can be obtained as T^(^ppg^ = n/ (r/j -f- J12/2) 
for j = 1, 2, 77, = 0, 1, 2, ... , where ijj = SjztJoj/2. We set 

61 ^ 62 ^ 400 in the experiment. For the entangled state, 

2 2 
Hm =27Tf2 Sjll + 27rJ'i2li -1^+271 J2 Jojl'^H, where 

2 

Hs = 2TrJ2 + 27rJi2li • I^ which is Heisenberg- 
i=i 

XXX-coupling (or isotropic Heisenberg) type. Similarly, 
'''m(ENT) = "/('^j >7oj/2). It is easy to find that |0+) 
is a non-degenerate eigenstate of i?s, and [Hg, Hint] = 
when 5i = 62 and Joi = Jm- In experiment, 5i = 62 = 



4 



(a) 




20 40 60 80 100 20 40 60 80 100 

Number of pulses Number of pulses 




FIG. 3: (color online). Experimental QZE with 
entanglement-based measurement on single qubit system, (a) 
Experimental observations of QZE with a series of RF pulses 
and measurements (left side) and with GRAPE engineer- 
ing (right side). The amplitude decay in the left pannel 
is simulated by adding an exponential decay exp( — fer™) 
to each Af-process (fc = 1.25, 0.8 and 0.7, respectively, 
for Tm ~ 2.517ms, 2.55ms and 3ms). (b) Evolution 
paths of the system on the Bloch sphere, corresponding to 
the three experiments with the measurement times Tm = 
2.517ms, 2.550ms, 3.000ms. 




FIG. 4: (color online). (a) Experimental QZE with 
entanglement-based measurement on two-qubit system for a 
product state (top plots) and (b) an entangled state (bottom 
plots), (b) Real part of the reconstructed density matrix of 
the initial and final state for preserving the entangled state 
via QZE. 



IV. FURTHER DISCUSSION 



200, Ji2 = 100, and Joi = J02 = 250. Note that J,j arc 
not the same as ones in the nature Hamiltonian Hnmr- 
Consequently, we can experimentally achieve the evolu- 
tion of Hm by quantum simulation technique [14] . 

However, the direct simulation of the Hamiltonian Hm 
requires a long operation time to realize the measure- 
ment M (T„i), especially in the case of the entangled state 
(about 26 ^ 31ms) (Fig. 2). Relaxation will be a seri- 
ous problem. Accordingly, we adopted the GRAPE en- 
gineering here for precise quantum control. The experi- 
mental results are shown in Fig. 4, which illustrates the 
Tm-dependent behavior in entanglement-based QZE both 
for the product state and the entangled state, like the 
single-qubit system. At the critical measurement time 
Tm e.g., r^(pps) = 3.059ms and t;:^(^ent) = 3.077ms, 
2? oscillates almost from to 1. When t,„ is right 
around r,* , e.g., t^{pps) = 3.1ms and t^(ent) = 3.2tos, 
we can see the amplitude of oscillations of V decays 
quickly. T) has almost no change at Tjn{pps) = 3.5tos, 
Tm{ENT) — 3.7ms, representing the QZE occurs and the 
preservation of quantum states works well. To assess the 
quality of the QZE, we also carried out full quantum state 
tomography for the initial state and the final state after 
N pulses for the case of the entangled state, shown in 
Fig. 4(c). With the help of the GRAPE pulses, we have 
successfully preserved the entangled state with a high fi- 
delity running up to 0.99 [15]. 



In the end, we further analyze the dynamic theory of 
QZE to explain the experimental results. For the single- 
qubit system, the first free evolution of U (r) evolves 
the initial state |s) = |0) to a superposition state: 
U{t)\0) = a\Q) + b\l). Then, it undergoes a mea- 
surement process M{Tm) to \ipi) — M{Tm)U (r) |0) = 

g-»7r(5i + ,/oimo)r,„^|0) + g»7r(5i-H Joi™o )r„ ^| ^j^j^ ^ 

±1/2 depending on a state of spin /q, |0) or |1). The 
measurement adds a phase shift between |0) or |1) of the 
system state: A0 = 27r((5i -I- Joimo)''"™. Therefore, we 
can see three different situations from the phase shift, 
(i) When = T^ = n/{5i + Joimo) (n = 0, 1, 2, . . . ), 
the phase shift A<?!) = 2n7r and = e"™'' (a|0) 4- &|1)). 
This implies that the measurement with the critical time 
only gives a whole phase factor e~*"'^ to the state. The 
repeated applications of the pulse and measurement drive 
the system to undergo a Rabi oscillation like the be- 
havior of a continuous RF pulse except for a whole 
phase factor, and no QZE occurs, [ii) When Tm ~ 
[n — l/2)/((5i + Joimo), the phase shift A0 = (2n — l)7r 
and iV'i) = e-''("-i/2)7r (•^|Q^) „ g^^^^^ acquires a 

phase of —1 with respect to |0) (that is, a 7r-phase shift). 
Hence, the evolution of the system spin is reversed after 
each measurement is applied, thus reversing the free evo- 
lution of the system and QZE occurs. This has a similar 
behavior to "bang-bang control" [9, 16]. It is not neces- 
sary to apply perfect 7r-phase shifts to lock the spin, (iii) 
When the measurement time has a small deviation ^ from 
T*-^, i.e., Tm = + ^, the intermediate case occurs. We 
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have = e-i27r(5i+Joimo)/]r™g-i27r(5i7i+Pi4)T|Q^) ^ 

g-i2n7r/ig-i27r{[(5i + Joimo)C/T+5i]/i-j27rPi/i}r|Q^ rpj^j^ 

functions as a RF rotation around an axis in the XZ 
plane, and the situation is similar to the first one, only 
with different amplitude. In fact, this approximation al- 
ways makes the whole unitary propagator [M{Tm)U{T)]^ 
for N repeated cycles good in the wide range of ^ in 
our experiment due to t — >■ 0. However, when \{Si + 
JqiTOo)?! ^ PiT, M(t„i)U{t) is a rotation almost around 
the z axis, which results in the initial state |0) unchanged. 
When 1^1 reaches 0.5/|(5i + Joimo\, we exactly return to 
the situation (ii). Consequently, the intermediate phe- 
nomenon occurs only in a range of \{Si + Joimo)^\ ^ Pit, 
where Pit « 0.018 in the experiment. Figure 3(b) shows 
clearly the corresponding evolutions of the system qubit 
on the Bloch sphere for these three situations. The sim- 
ilar analysis can be also used to explain the experiments 
on the two-qubit composite system. 



and the composite system. In the single-qubit system, 
the free evolution Hamiltonian is Hfree = 27r((5i/^-|-Pi/^) 
and the measurement Hamiltonian = 2t:{5iII + 

Joil^ll). Following the same procedure in Ref.[7], the 
total time evolution operator is determined as 



Utot {t) 



exp 



'2TTiSillt 



Pit 



M (NTm) . 



Here, we neglect the high order terms of r in the exponent 
and denote that 



+ sin[7rr„(^i + Joi/O)] 



V. CONCLUSIONS 

In conclusion, we have experimentally demonstrated 
the QZE with entanglement-based measurement in a 
NMR ensemble, where the frequent measurements are 
implemented through the dispersive-coupling between 
the target system and the measuring qubit. Therefore, 
the measurement is dynamically described by the uni- 
tary evolution of duration Tm, rather than the projec- 
tion measurements. With a three-qubit NMR system, 
we have successfully observed the dynamic QZE of a 
single-qubit system as well as the composite system, espe- 
cially an entanglement preservation. Our experiment also 
clearly shows the dependence of this dynamic QZE on 
the measurement time in the frequent measurements: 
the system exhibits no QZE at certain critical measure- 
ment times T^. This well distinguishes from the usual 
one based on projection measurements. Moreover, we 
have experimentally demonstrated nontrivial quantum- 
state steering towards the efficient preservation of entan- 
glement using the dynamical QZE. 



where l\_ = ll 



ily. 



When the measurement time 
approaches to the critical value t^ = n/ {5i + Joimo), 
[n = 0,1,2,..., and mo = -1-1/2 or —1/2, the eigen- 
value of /° according to its the initial state) , /i is propor- 
tional to Nil, tlius Utot {t) ~ exp [-27ri(5i/i + Pill)t] 
will drive the system oscillating and the QZE is violated. 
Otherwise, when the measurement time is not around 
T*j, /i is finite hence Utot (t) becomes a unit operator in 
the limit r — >■ 0. 

For the product state |00) of the composite system, 
the measurement Hamiltonian with Ising-type coupling 

2 2 

is Hm ^ 2nJ2 +2tt J2 JjkPjt We find that 

i=i i<fc,=o 
the total time evolution operator reads 



Utot {t) « 
exp|-iEj=i,2 

where 



2TT6jIit + ^ [fj + h.. 



} M [Nt^) , 
(Al) 
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Appendix A: critical measurement times 

Now we analyze the critical measurement times for two 
typical cases in the experiments: the single-qubit system 



-J ^-»^(jV+l)r^(r,, + lJi.) Sin [irNTm (% + \Jl2)] 

+ sin [-KTyn {V] + ^^12)] 



where = + il^ and rjj = 5j + JojTOo- When 
the measurement time approaches to the critical value 
T*j = n/ {f]j + J12/2), (n ~ 0,1,2,...), as the same 
reason in the single-qubit case, Utot (t) is proportional 
to exp[—27ri{Sjli+Pjll)tj, which violates the QZE. 
Similarly, we can obtain the critical measurement time 
T*j = 71/ {5j ± Jaj/2) for the case of the XXX coupling. 
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